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The terahertz (THz) frequency range (0.1-10 THz) fills the gap between the microwave 
and optical parts of the electromagnetic spectrum. Recent progress in the generation 

i n 

(UD and detection yD of the THz radiation has made it apowerful tool for fundamen- 
tal research and resulted in a number of applications J^; 0). However, some important 
components necessary to effectively manipulate THz radiation are still missing. In par- 
ticular, active polarization and phase control over a broad THz band would have major 
applications in science and technology. It would, e.g., enable high-speed modulation for 
wireless communications J3) and real-time chiral structure spectroscopy of proteins and 
DNA dfjl; ItJ). In physics, this technology can be also used to precisely measure very weak 
Faraday and Kerr effects, as required, for instance, to probe the electrodynamics of topo- 
logical insulators dg; |9|). Phase control of THz radiation has been demonstrated using 
various approaches. They depend either on thephysical dimensions of the phase plate 



(and hence provide a fixed phase shift) (|10i : 
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121) or on a mechanically controlled time 
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delay between optical pulses (and hence prevent fast modulation) (1131 : 114c Il5ft . Here, we 
present data that demonstrate the room temperature giant Faraday effect in HgTe (|16|) 
can be electrically tuned over a wide frequency range (0.1-1 THz). The principle of op- 
eration is based on the field effect in a thin HgTe semimetal film. These findings together 
with the low scattering rate in HgTe open a new approach for high-speed amplitude and 
phase modulation in the THz frequency range. 

The idea behind our experiments is to exploit the strong interaction between a solid-state 
plasma and THz radiation (|17r) . Because of this interaction, left- and right-handed circularly 
pola rized THz waves have different complex refraction indices when a magnetic field is applied 



(|18r) . Recently, this characteristic property has been proposed to realize broadband THz mod 



ulators and phase plates based on electron-doped InSb crystals (19). However, there are two 
problems in the realization of such devices for practical applications. First, the giant Faraday 
effect in InSb does not survive at room temperature. Second, this approach requires fast modu- 
lation of the applied magnetic field with an amplitude of several hundreds mT, which is a rather 
challenging demand. 

In semimetals with zero (or very small) band gap, such as HgTe, the electron plasma exists 
in undoped samples, due to the thermal activation of electrons. The ability to use such undoped 
HgTe layers results in a very high electron mobility and, consequently, we have been able to 



demonstrate the giant Faraday effect even at room temperature (|16l) . 

To achieve control over the Faraday effect in a constant magnetic field, we have now fabri- 
cated devices fitted with transparent gate electrodes, which allow to change the carrier density 
and thus the properties of the electron plasma in the HgTe layer. By applying a moderate volt- 
age \U\ < 1 V to these gates we are able to control the Faraday rotation and Faraday ellipticity, 
i.e., amplitude, phase shift, and polarization of the THz radiation at room temperature. 

The THz experiments have been carried out in Faraday geometry as shown schematically 
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Figure 1 : Electric voltage control of THz radiation, (a) - Schematic representation of exper- 
imental arrangement to control the Faraday rotation (9) and Faraday ellipticity (77) by a gate 
voltage. (b,c) - Normalized detector signal as function of gate voltage in a geometry with par- 
allel (b) and crossed (c) polarizers, respectively. 
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in Fig. la. In this geometry the external magnetic field is parallel to the propagation direction 
of the THz radiation. The incident polarization is linear and the polarization of the transmit- 
ted radiation is detected by a wire grid analyzer. Further experimental details are given in the 
Supplementary Material. All experiments presented here have been carried out at room tem- 
perature. Our estimations suggest that the concentration of carriers confined to the topological 



201) is comparable with the concentration of carriers occupied 



surface state in this material 
bulk states. 

As shown in Fig. 1 be, we observe a substantial variation of the amplitude of the transmitted 
signal as a function of gate voltage, both in parallel (t p ) and in crossed (t c ) polarizer geometry. 
Consequently, the Faraday rotation and ellipticity can be controlled by using a moderate electric 
voltage. In a first approximation, the observed changes may be attributed to the gate dependence 
of the effective charge density in the HgTe film. As will be shown below, not only the density, 
but all electrodynamic parameters of the HgTe are affected by the gate voltage. We mention 
already at this point that the response of the THz signal is asymmetric with the sign of the 
gate voltage, which is probably due to carrier trapping at the interface between HgTe and the 
gate insulator. We also observed hysteresis effect when tuning the gate voltage beyond certain 
critical value. 

In order to characterize the electron dynamics in HgTe under applied gate voltage, a set of 
experiments at different magnetic fields B and THz frequencies cu has been carried out. Figure 
2 demonstrates the magnetic field dependencies of the magnetooptical signal of the HgTe film at 
three different gate voltages. The transmission amplitude for parallel polarizers (Fig. 2a) shows 
a clear minimum at a frequency-dependent field position which corresponds to the cyclotron 
resonance in HgTe f2 c = to. The width of the minimum in t p is directly connected to the 
scattering rate of electrons. As the width in Fig. 2a changes strongly with gate voltage, a 
variation of the characteristic scattering time r is evident. 
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Figure 2: Magnetooptic transmittance of the HgTe film. Magnetic field dependent amplitude of 
the transmitted radiation in Faraday geometry and for different gate voltages as indicated, (a) 
- Transmission in parallel polarizers geometry. The arrow marks the position of the cyclotron 
resonance cu = f2 c . (b) - Transmission in crossed polarizers. The maximum indicated by the 
arrow allows a determination of the quasiparticle scattering time r. 
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Figure 2b shows the amplitude of the transmission in crossed polarizer geometry. All curves 
reveal a characteristic maximum at a magnetic field around ±2 T. As shown in the Supplemen- 
tary material [Eq. (fTOh l. this maximum corresponds to a field value of VL c t = 1 and it allows 
for an independent determination of r. Again, a strong variation of the maximum is observed in 
the field dependence of t c , which confirms the variation of the scattering time with gate voltage. 

In addition to the transmission amplitude, phase information has been obtained using a 
Mach-Zehnder interferometric arrangement. From these experiments, the polarization state of 
the transmitted THz radiation can be fully polarized, including Faraday rotation and ellipticity. 
These data are shown in Fig. 3 for the parameter range accessible in our setup (\B\ < 7T and 
\U\ < 0.8 V). The data in Fig. 3 are given for two different frequencies, close to the upper and 
the lower limits of our spectrometer. From the Faraday ellipticity vs. Faraday rotation diagram 
of Fig. 3 one can see that a desired polarization state can be set by varying either magnetic field 
or gate voltage. 

The observed gate dependence of the magnetooptical spectra can be well interpreted within 
the Drude model. The details of the characterization procedure and extraction of the electro- 
dynamic parameters are given in the Supplementary Material. In brief, we fit simultaneously 
the magnetic field dependencies using charge density (n 2D ), scattering rate (1/r) and cyclotron 
frequency (in terms of Vt c /B) as the only fitting parameters. By repeating this fit for different 
gate voltages, we determine how these parameters depend on U, as shown in Fig. 4. Two- 
dimensional electron conductivity derived from the other parameters as 

a 2D = en 2D —r (1) 

is also shown in Fig. 4. 

As expected, a strong variation of the electron density with gate voltage (nearly a factor of 
two) is observed (Fig. 4a). Assuming the simple capacitor formula An 2D = £o^Si 3 N 4 A£/ydsi 3 N 4 
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Faraday Rotation (°) 

Figure 3: Tunability of Faraday rotation and ellipticity. Parametric presentation of the attainable 
range of Faraday rotation and ellipticity for different fields and voltages at 0.98 THz (a) and 
0.18 THz (b), respectively. The data are shown as magnetic field dependent polarization states 
at three different gate voltages, as indicated. The values of the external magnetic field are 
color-coded. Symbols - experiment, lines are a guide to the eye. 
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Figure 4: Gate dependence of the electrodynamic parameters. Electrodynamic parameters of 
the conduction electrons in HgTe at room temperature as obtained from Faraday rotation and 
ellipticity. a - Effective two-dimensional electron density, b - Two-dimensional dc conductivity, 
c - Electron scattering rate, d - Cyclotron resonance frequency. 
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to estimate the density variation, we obtain An 2 D — 4.15 • 10 12 cm -2 per Volt. Here 6q is the 
permittivity of the free space and £si 3 N 4 = 7.5 is the dielectric constant of SiaN 4 insulating 
layer. This is in a reasonable agreement with our experimentally observed value of 1.12 ■ 10 12 
cm~ 2 V _1 . Where again the difference with the capacitor model likely results from carrier trap- 
ping at the HgTe/Si 3 N 4 interface. 

Contrary to the increase of the electron density with the gate voltage, the conductivity (ct 2 _d, 
Fig. 4b) is a decreasing function of U. This is rather unusual and can conceivably be explained 
by increased carrier scattering (r _1 , Fig. 4c) at the HgTe/Si3N 4 interface for larger gate volt- 
ages. Finally, we observe that the cyclotron frequency depends on the gate voltage as well, Fig. 
4d. This result is expected for the Dirac surface states in HgTe with Q, C J B oc 1/ v /ra 2 D- It can 
be also caused by the deviation from ideal parabolic dispersion of bulk electrons in HgTe. 

All results presented here have been obtained at room temperature, suggesting practical 
applications, such as fast phase and amplitude modulation or direct control of the polarization 
state by gate voltage and/or magnetic field. As the observed tuning of the THz radiation is 
due to modulation of the electron conductivity in the HgTe layer and is obtained by purely 
electrical means, it should be possible to achieve rather high modulation speed, comparable to 
this obtained with high electron mobility transistors (HEMT). The modulation amplitude is as 
large as several degrees per Volt, and we believe that it can be further improved by optimizing 
parameters of the gate material, insulating barrier and HgTe layer. 
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Supporting Online Material 

The sample studied in this work is a coherently strained 100-nm-thick nominally undoped HgTe 
layer, grown by molecular beam epitaxy on an insulating CdTe substrate (|2ll). A 3 nm thin 
R11O2 layer has been used as a gate electrode and 10 nm thick layer as an insulating 

barrier. The gate electrode transmits approximately 80 % of the terahertz radiation and revealed 
no own Faraday rotation signal. The amplitude of the gate voltage has been limited to ±0.8 V 
to avoid an electrical breakdown of the insulating layer. 



Transmittance experiments at terahertz frequencies (0.1 



iz < v < 1 THz) have been 



23|) which allows measurement 



carried out in a Mach-Zehnder interferometer arrangement (|22l : 
of the amplitude and phase shift of the electromagnetic radiation in a geometry with controlled 
polarization. Using wire grid polarizers, the complex transmission coefficient can be obtained 
both in parallel and crossed polarizers geometry. Static magnetic fields, up to ±7 Tesla, have 
been applied to the sample using a split-coil superconducting magnet. 

To interpret the experimental data we use the ac conductivity tensor obtained in the 
classical (Drude) limit from the Kubo conductivity of topological surface states (see e.g. Ref. 



(|24[) ). The diagonal, a xx {uS), and Hall, a xy {u), components of the conductivity tensor as func- 



tions of THz frequency uj can be written as: 



1 — tUJT 



a x M = -a yx (u) = {1 _ lUT)2 + {QcT)2 *o ■ (3) 

Here f2 c = eBv F /hkp is the cyclotron frequency, (Tq is the dc conductivity, B is the magnetic 
field, Vp, hp, e, and r are the Fermi velocity, Fermi wave-number, charge, and scattering time 
of the carriers, respectively. For the Dirac spin-helical surface states the Fermi wave-number 



depends on the 2D carrier density, n 2 o, through relation k F = \j^Kn%r>, with no spin degen- 
eracy. For massive carriers the cyclotron resonance frequency can be written as f] c = eB/m, 
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where m is the effective electron mass in the parabolic approximation. 

For a free standing film and neglecting any substrate effects, the complex transmission co- 
efficients in parallel (t p ) and crossed (t c ) polarizers geometry can be written as: 

t = 4 + 2S - (4 ) 

tp 4 + AE XX + E 2 XX + El y ' {) 

t = '—^L (5) 

c 4 + AE XX + E 2 XX + ' 

Here E xx and E xy are effective dimensionless 2D conductivities, defined as: E xx = a xx dZ 

and E xy = a xy dZ with the HgTe film thickness d = 100 nm and the vacuum impedance 



Z w 377 tt. 



formalism ([25 



n order to take into account the influence of the substrate a transfer matrix 



26; 



161) is utilized. The electrodynamic properties of the CdTe substrate are 
obtained in a separate experiment on a bare substrate. 

According to Eq. ©, the transmittance in parallel polarizers (t p ) depends mainly on E xx . 
The minima in \t p \ roughly correspond to the cyclotron resonance energy and they scale with 
magnetic field. This may be understood by taking into account that for our film E 1 and Eqs. 
(|4I5T) simplify to: 

tp — 1 — E xx /2] t c ~ E xy /2 . (6) 

In the limit cut 3> 1, Eq. © may be approximated by 

1 - iur 

which leads to a resonance like feature in a xx and t p for fi c = oj. Thus, the positions and widths 
of the minima in Fig. 2a are directly connected to Fermi velocity (vp) for linear dispersion [or 
the effective mass (m) for parabolic dispersion] and the scattering rate (r" 1 ) of the electrons. 

The polarization rotation 9 and the ellipticity rj are obtained from the transmission data 
using: 

tan(20) = 2»(x)/(l - |x| 2 ) , (8) 
11 



sm(2 V ) = 29f(x)/(l + |x| 2 ) • (9) 

Here x — t c /t p and the definitions of 9 + ir] are shown schematically in Fig. la. A direct 
interpretation of the complex Faraday angle is in general not possible because of the interplay 
of a xx and a xy in the data. 

In the low frequency limit (lot < 1) Eq. © simplifies to the static result: 

**» = i+%ry a ° ■ m 

The last expression has a maximum at Vl c (B)t = 1, which leads to maxima in t c (Eqs. ([5] [6]) 
and Fig. 2b) and in 9 at about the same field value. Therefore, the Faraday angle provides a 
direct and an independent way of obtaining the scattering rate 1/r. 
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